Weak Lensing of Galaxy Clusters in MOND 
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ABSTRACT 

We study weak gravitational lensing of galaxy clusters in terms of the MOND (MOdified 
Newtonian Dynamics) theory. We calculate shears and convergences of background galaxies for 
three clusters (A1689, CL0024+1654, CL1358+6245) and the mean profile of 42 SDSS (Sloan 
Digital Sky Survey) clusters and compare them with observational data. The mass profile is 
modeled as a sum of X-ray gas, galaxies and dark halo. For the shear as a function of the angular 
radius, MOND predicts a shallower slope than the data irrespective of the critical acceleration 
parameter go. The dark halo is necessary to explain the data for any go and for three interpolation 
functions. If the dark halo is composed of massive neutrinos, its mass should be heavier than 
2 eV. However the constraint still depends on the dark halo model and there are systematic 
uncertainties, and hence the more careful study is necessary to put a stringent constraint. 



Subject headings: cosmology: theory - 
CL1358+6245) - gravitation - gravitational lensing 

1. Introduction 

MOND (MOdified Newtonian Dynamics) Q is 
a theoretical alternative to Newtonian dynamics, 
proposed by Milgrom (1983). The theory itself 
strengthens gravitational force at large distances 
(or small accelerations) to explain galactic dynam- 
ics without dark matter. The equation of motion is 
changed if the acceleration is lower than the criti- 
cal value go ~ 1 x 10~ 8 cm/s 2 . It is well known that 
this theory can explain galactic rotation curves 
with only one free parameter: the mass-to-light ra- 
tio (see review Sanders & McGaugh 2002). There 
are two motivations to study such an alternative 
theory: (i) General Relativity (GR) has not been 
tested accurately at much larger scale than 1 AU 
(ii) dark matter particles have not been directly 
detected and their nature still eludes us. Under 
these circumstances, several authors have recently 
studied alternative theories to GR (e.g. Aguirre 
2003). 

Bekenstein (2004) recently proposed a rel- 



axies: clusters: individual (A1689, CL0024+1654, 



ativistic covariant formula of MOND (called 
TeVeS) by introducing several new fields and 
parameters. Following this, many authors be- 
gan discussing relativistic phenomena such as 
parameterised Post-Newtonian f ormalism in the 
solar system (jBekenstein 20041 ). gravitational 
lensing (e.g. Zhao et al. 2006), cosmic mi- 
crowave backgr ound and large scale structure 
of the Universe dSlosar et al. 20051: ISkordis 20051: 
Skordis et al. 20061: iDodelson fc Ligu ori 2006!). In 



lr The phrase MOND is used here to refer to modified New- 
tonian gravity models without any dark matter. 



this paper, we discuss weak gravitational lensing 
of galaxy clusters. 

Weak lensing provides an important observa- 
tional method with which to test MOND. This 
is because weak lensing probes the lens poten- 
tial outside of the Einstein radius, te ~ (AID) 1 / 2 
~ 150kpc(Af/10 14 M Q ) 1 / 2 (L>/ J ff,7 1 ) 1 / 2 , where M 
is the lens mass and D is the distance to the source. 
The gravitational law changes outside the MOND 
radius, ru = (M/go) 1 / 2 , where the acceleration is 
less than go. Since go ~ Ho/6, the Einstein radius 

is a few times smaller than the MOND radius 
r M, by a factor 2 at least, in the cosmological sit- 
uation. Hence, we can test the MOND-gravity 
regime by weak lensing. 
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Weak lensing is superior to X-rays as a means 
of probing the outer region of clusters. The lensing 
signal (strength of the shear) is proportional to the 
surface density E. On the other hand, the X-ray 
luminosity is proportional to the density squared 
p 2 and hence X-rays can probe the inner regions 
of clusters. Hence, the outer region of the clusters 
can be probed with weak lensing. 

The gravitational lensing in MOND has been 
studied by many authors. Before Bekenstein 
proposed the relativistic formula, some assump- 
tions were madjl to c alculate the lensing quan- 
tities (iQin et al. 19951: iMortlock fe Turner 2001a: 



2001bt White fe Kochanek 200lL iGavazzi 20021 ). 
Just after Bekenstein's proposal, Chiu et al. 

(2006) and Zhao et al. (2006) first studied the 
lensing in detail and tested MOND with strong 
lensing data of galaxies. Zhao and his col- 
labo rators studied the gravitational lens statis- 
tics ( Chen fe Zhao 2006h and investigated a non- 
spherical symmetric lens (jAngus et al. 20061 ) . Re- 
cently, Clowe et al. (2006) indicated that a merg- 
ing cluster IE 0657-558 cannot be explained by 
MOND because the weak lensing mass peak is 8a 
spatial offset from the baryonic peak (= mass peak 
of X-ray gas). However, Angus et al. (2007) noted 
that MOND can explain the data if the neutrino 
halo is included (see also Feix et al. 2007). Fur- 
thermore, the high observed collision velocity of 
the bullet clusters (shock velocity of ~ 4700km/s) 
i s more readily obtained in MOND than CDM 
(I Angus fc McGaueh 20071) . 

Jee et al. (2007) recently found a ring like dark 
matter structure at 9 ~ 75" in CL0024+1654 by 
analysing strong and weak lensing data. They sug- 
gested that the ring was formed by the line-of-sight 
collision of two clusters, like the bullet cluster IE 
0657-558. Just after their discovery, Famaey et al. 

(2007) noted that MOND can easily explain the 
density of the ring by adding the massive neutrino 
of 2 eV. 

In this paper, we study three clusters (A1689, 
CL0024+1654, CL1358+6245) and mean profile 
of 42 SDSS (Sloan Digital Sky Survey) clusters. 
We calculate shears and convergences for these 
clusters and compare them with the observational 



2 For example, Qin et al. (1995) assumed that the bending 
angle is 2 times larger than that for massive particles in the 
limit of m — > by analogy with GR. 



data. We perform a \ 2 fit of the data to give a 
constraint on the dark halo profile and the neu- 
trino mass. Throughout this paper, we use the 
units of c = G = 1. 

2. Basics 

We briefly review the basics of gravitational 
lensing based on the relativistic MOND theory for 
a spherically symmetric lens model. Detailed dis- 
cussions are given in Bekenstein (2004) and Zhao 
et al. (2006). 

When a light ray passes through a lens with the 
impact parameter b, the deflection angle is 



a(b) = 2b 



dl 



sO) 



(1) 



where I is the distance along the light path and r 
is the distance from the lens center, r = \Jl 2 + b 2 
(see FigO}. 

The gravitational force due to the lens is 

Hg/9o)g(r)=g N (r) = ^^, (2) 

where is the usual Newtonian acceleration and 
M{< r) is the lens mass enclosed inside a ra- 
dius r. We use a standard interpolation function^] 
jj,(x) = xj\J\ + x 2 with go = 1 x 10~ 8 cm s~ 2 . 
Then, jl(x) = 1 (i.e. g = <7at) for g 3> 50, while 
p{x) = x (i.e. g = ^fg~NgV) for g < g . 
The lens equation is 



D 



LS 



a(8). 



(3) 



Here, 8 S and 6(— b/ Dl) are the angular positions 
of the source and the image, and Dl, D$ and D^s 



3 We also examine other interpolation functions to study its 
dependence in section 4. 

light path / 



lens 




Fig. 1. — A schematic picture of the light ray pass- 
ing through the lens. 
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are the angular diameter distances between the 
observer, the lens and the source@ The shear 7 
and the convergence k are given by, 



7= 2 



d9 



K =2 



d6 



(4) 



We note that if the mass increases as M oc r p 
withp > 0, the shear and the convergence decrease 
as 



7 oc k oc 9 P 2 
oc OP/ 2 " 1 



for g > 50, 
for g < g , 



(5) 



from Eqs.|T|)-(I3]). The slopes of 7 and k for g ^> go 
are steeper than that for 3 <C go- This is because 
the gravitational force is proportional to g 1 ^ 2 for 
g <C go, an( l hence the force decreases more slowly 
at larger distances. Comparing the slope in Eq.© 
with the observational data, we can test MOND. 

3. Analysis with Cluster Data 

We calculate the shear 7 and the convergence k 
based on the MOND theory for the three clus- 
ters, A1689, CL0024+1654, CL1358+6245, and 
the mean profile of 42 SDSS clusters. The mass 
profiles of these clusters have been measured by 
gravitational lensing for a wider range of angular 
diameters, and hence these clusters are an appro- 
priate system to investigate the angle-dependence 
of the shear and the convergence. 

3.1. A1689 

Several authors have been studying the mass 
profile of the rich cluster A1689 at 2 = 0.183^ by 
strong and weak lensing, X-ray emission of gas, 
and dynamics of cluster members (e.g. Limousin 
et al. 2006 and references therein). The analysis of 
lensing data shows a small ellipticity (e = 0.06 in 
Halkola et al. 2006) and supports the assumption 
of quasi-circular symmetry (Umetsu, Takada & 
Broadhurst 2007). Andersson & Madejski (2004) 
provided the hot gas mass profile (40kpc < r < 
IMpc) directly determined by X-ray observational 



4 We use the distance in the usual FRW (Friedmann- 
Robertson- Walker) model with Qm = 1 — = 0.3 and 
H = 70km/s/Mpc. The distance in MOND is alm ost the 
same as that in the FRW model llBekenstein 2004T) . 

5 l' corresponds to 184kpc. 



data of the XMM-Newton telescope. Zekser et al. 
(2006) gave the galaxy mass profile (20kpc < r < 
260kpc) from the surface brightness profile, as- 
suming the constant mass-to-light ratio 8M Q /L© 
(B-band). FigJ^a) shows the mass profiles of the 
gas (dotted), the galaxies (dashed), and the sum 
of them (solid line). We also show the dark halo 
profile which will be needed to match the observa- 
tional data (we will discuss this later). 

FigEKb) shows the Newtonian gravitational ac- 
celeration gjv normalized to go- As shown in 
this panel, the transition radius corresponding to 
9N = 9a (denoted by a horizontal dotted line) is 
100 kpc for the gas + galaxies and is at 1000 kpc 
if the dark halo is added. 

Broadhurst et al. (2005a) measured the distor- 
tions of 6000 red galaxies over 1' < 9 < 20' by the 
wide field camera, Suprime-Cam, of the Subaru 
telescope. Panel (c) shows their results, reduced 
shear profile 7/(1— k) ■ The mean source redshift is 
z s — 1±0.1 based on a photo-z estimation for deep 
field data. The solid line is the MOND theoretical 
prediction with z a = 1. The gravitational source 
is only baryonic component (gas + galaxies). We 
note that for 9 < 10' the solid line is clearly smaller 
than the data. This indicates that the gravita- 
tional force is too weak to explain the data. In 
order to solve this discrepancy, we need a very 
high mass-to-right ratio ~ 200Mq/L^. Even if 
the critical acceleration value go increases, the dis- 
crepancy could not be resolved. In this case, the 
amplitude of the shear increases but the slope is 
too shallow to fit the data. MOND predicts shal- 
lower slope than 7 oc 9~ 1 for g < go (since p > in 
Eq. (JSJ) ) , while the data in panel (c) clearly shows a 
steeper slope than this. Hence MOND cannot ex- 
plain the data for any mass model and any acceler- 
ation parameter go in the low acceleration region 
9 < So- 

We comment on the dependence of the source 
redshift z s on the above results. The quantities 
7 and n depend on z s through a combination of 
Dls/Ds, 7 oc k oc D LS /D S , from Eqs.© and 
((4]). Hence these slopes are independent of z s and 
the above results do not change. Furthermore the 
quantity Dls/Ds is not sensitive to z s for rela- 



° The shear and the convergence are proportional to the 
mass-to-light ratio M/L for g go, while (M/L) 1 / 2 for 
9 < 90- 
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tively low lens redshift (z — 0.183 for this cluster). 

Previously, Aguirre, Schaye & Quataert (2001), 
Sanders (2003), and Pointecouteau & Silk (2005) 
reached the same conclusion as ours by studying 
temperature profiles of clusters. They indicated 
that the temperature data near the core is higher 
than the MOND prediction. Sanders (2003) noted 
that if the dark matter core were added, this dis- 
crepancy could be resolved. Following the previ- 
ous studies, we include the dark halo to explain 
the observational data. We use the dark halo with 
a flat core : 



M(< r) = M 



'•o 



(0) 



Here r is the core radius and the density steeply 



decreases with proportional to 



for 



r > r Q 



: p oc (r + r )~ 4 . We perform a x 2 nt 01 the 
data in order to determine the parameters Mq 
and r . We also use the strong lens data (conver- 
gence field k) from Broadhurst (2005b) in panel 
(d) in order to put a strong constraint on the 
core radius. They constructed the convergence 
field at 0.08' < 9 < 1.4' from 106 multiply im- 
ages of 30 background galaxies by Hubble Space 
Telescope Advanced Camera for Surveys. The x 2 



is given by x = J2 t ( x i 



~data\2 



) 2 /°f where Xi is 



the reduced shear 7/(1 — k) and the convergence 
k at the i-th angle, xf ata is the data and oi is 
the standard deviation. The best fitted model is 
M = (l.l±0.1)x 1O 15 M and r Q = 174±11 kpc. 
These relative errors are less than 10% because of 
combining the strong and weak lensing data. The 
minimized % 2 -value per degree of freedom (dof) 
is Xmin/dof = 14.7/26. The results are insensi- 
tive to the mass-to-light ratio. As shown in panels 
(c) and (d), this model (dashed line) fits the data 
well. The dashed line in panel (c) is steeper than 
the solid line, since 9 < 5' («-> r < lOOOkpc) is 
the high acceleration region g > go from panel (b) 
and hence the slope is steeper, as can be seen in 
Eq.©. 

We also try to fit the data by the other halo 
profiles : Navarro-Frenk- White (NFW) model p oc 
r(r + ?"o)~ 2 7 Hernquist model p oc r(r + ro)~ 3 , 
isothermal (IS) with a core p oc l/(r 2 + r 2 ,). The 
minimum Xmin are 29-8 for NFW, 22.6 for Hern- 
quist, and 49.6 for IS with core which are larger 
than 14.7 for our model in Eq.©. This is be- 
cause the convergence data favor a flat core and 



the shear data favor a steeply decreasing profile. 

3.2. CL0024+1654 

CL0024 is a rich cluster at z = 0.395 (1' cor- 
responds to 32 kpc) with a velocity dispersion 
of 1200km s^ 1 (iDressler fc Gunn 1992ft . FigH (a) 
shows the mass profiles of the gas determined b y 
the XMM-Newton telescope (jZhang et al. 20051 ). 
the galaxi es with a mass-to- right ratio 8M /L 
(K-band) (jKneib et al. 2003 ). and the dark halo. 
For the larger radius > 2 Mpc the baryonic mass 
exceeds the dark halo mass. This is because we 
extrapolate the gas profile (fitted by isothermal /3 
model for r < 1 Mpc) to the larger radius. Kneib 
et al. (2003) provided the reduced shear profile 
up to 10' measured by the Hubble Space Tele- 
scope as shown in panel (c). The mean source 
redshift is z s — 1.15. We also use a constraint 
from an angular position of Einstein radius at 27" 
(denoted by a black square ■) based on an obser- 
yation of giant arc s of distant galaxy at z = 1.675 
(|Smail et al. 19961 : iBroadhurst et al. 200ol ). The 
angular position of arcs is used to set a constraint 
on the enclosed mass within 27" (— 140 kpc). Sim- 
ilar to the previous case of A1689, for 9 < 10' the 
solid line is smaller than the data, and we need the 
~ 3OM0/L0 to solve this discrepancy. The best 
fitted halo model is M = (2.0 ± 0.4) x 10 14 M o 
and r = 42 ± 15 kpc in Eq.© with x^^/dof 
= 5.6/7. We caution that the core radius in the 
best fitting model ro = 42 ± 15 is smaller than the 
inner most data point (27" = 140 kpc) , and hence 
this result has little meaning. It only means that 
the enclosed mass is ~ Mq inside of 140 kpc. 

3.3. CL1358+6245 

The redshift of CL1358 is z = 0.33 and 1' cor- 
responds to 280 kpc. Fig Eta) shows the mas s dis- 
tribution of the gas ( Arabadjis et al. 20 02). the 
galaxie s with a mass-to-righ t ratio 8M Q /L (V- 
band) ( Hoekstra et al. 19981) . and the dark halo. 
Hoekstra et al. (1998) presented the reduced shear 
profile (9 < 4') measured by HST as shown in 
panel (c). They fitted the data by the isother- 
mal sphere model with the velocity dispersion of 
780 ± 50 km/s. The solid line is the MOND 
prediction with Dls/Ds — 0.62. We need the 
~ 30M Q /L Q to fit the data if we assume only 
baryonic components. The best fitted model is 
Mo = (6.9 ± 2.8) x 10 13 M Q and with r = 69 ± 31 
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(arcmin) (arcmin) 

Fig. 2. — Results for the cluster A1689. The top left panel (a): The mass profiles of the gas (dotted 
line), the galaxies (dot-dashed line), the gas + galaxies (solid line), and the dark halo (dashed line). The 
quantity M(< r) is the mass enclosed within the radius r. The top right panel (b): The Newtonian 
gravitational acceleration gw normalised to go for only baryonic components (gas + galaxies) (solid line) and 
all components (dark halo is added) (dashed line). The left bottom panel (c): The reduced shear 7/(1 — k) 
as a function of the angular radius. The data is from Broadhurst et al. (2005a). The solid line is the MOND 
prediction. On the dashed line, the dark halo is added. The right bottom panel (d): The convergence field 
k from Broadhurst (2005b). From panels (c) and (d), the MOND cannot explain the data unless the dark 
halo is added, because the gravitational force is too weak near the core. 
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Fig. 3. — Same as Fig[2j but for the cluster CL0024+1654. The black square ■ in panel (c) denotes the 
angular position of the Einstein ring. 
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kpc with Xmm/dof = 4.5/7. Same as the case for 
CL0024, the core radius is smaller than the inner 
most data point. Although the discrepancy be- 
tween the MOND prediction and the data is not 
so large in comparison with the previous cases, the 
dark halo model is better. 

3.4. SDSS clusters 

Sheldon et al. (2001) studied weak lensing of 42 
clusters in SDSS data. They provided the mean 
shear of 42 clusters up to the radius of 2700 kpc as 
shown in FigO The vertical axis is £ cr 7, where 
S cr = 1/ {4-k)Ds/ (DlDls) is the critical surface 
density, and the horizontal axis is the projected 
radiuqj- Here, S cr 7 does not depend on the source 
redshift. The data is well fitted by a powe r law 
with index -0.9 ± 0.3 ([Sheldon et al. 200l[ ). On 
the solid and dotted lines, we consider only the gas 
component, given by the isothermal beta model: 



PO) = Pa 



-3/3/2 



(7) 



with (3 = 0.6, r c = 100 kpc, p /p cr = 8000 (solid 
line) and 800 (dotted line), here p cr is the critical 
density at the present. As shown in the figure, the 
fit is poor. This is because the slope of the shear 
is -3/3/2 + 1/2 = -0.4 for g < g and it is flatter 
than the data. The dashed line is the dark halo 
model (given in Eq.([6])). The best fitted model is 
M = (3.8 ± 1.1) x 10 13 M o and r = 75 ± 34kpc 
with xL„/dof = 10.0/7. The dashed line fits the 
data well. 

4. Limit on Neutrino Mass 

In previous studies, several authors assumed a 
massive neutrino with a mass of ~ 2 eV as the 
dark matter to explain the observational data (e.g. 
Sanders 2003; Skordis et al. 2006). In this section, 
we put a constraint on its mass from the weak 
lensing of clusters. 

The neutrino oscillation experiments provide 
the mass differences between different species : 
AmJ < 10 _3 eV 2 . Here we consider massive neu- 
trinos whose masses are much heavier than Atoj, 



7 The quantity S C r7 is related to the surface density of the 
lens in GR : £ cr 7 = £(< R) - T,(R). 



and assume they are degenerate: they have (al- 
most) the same mass, independent of species. Us- 
ing the maximum phase space density h~ 3 , the 
maximum density of the neutrino dark halo is 
given by (Tremaine & Gunn 1979; Sanders 2003; 
2007) i 



= 2.3 x lO~ 5 M /pc : 



\2eVJ \keV J 



3/2 



(8) 

where T is the X-ray temperature : 9.0 ± 0.1 keV 
for A1689 (lAndersson fc Madeiski 200l . 3.5 ±0.2 
keV for CL0024 dZhang et al. 20051). an d 7.2 ±0.1 
keV for CL1358 (|Arabadiis et al. 20021 ). 

For A1689, the core density of the neutrino halo 
is p c = 3AV(47rr 3 ) from Eq.©. For CL0024 and 
CL1358, as we noted, the core radius ro in the 
best fitting model is smaller than the inner most 
data point. Hence, in order to put a conserva- 
tive bound, we use the mean density inside the 
second innermost data point r 2 , pi = p(< r%) 
= 3M(< r 2 )/(47rr|). Since p2 < p c , we obtain 
a lower bound of m„ from Eq.([8]). The results of 
p c and p2 are shown in Table 1. The lower row 
in each cluster is the case of another halo model 
M(< r) = r 3 /(r 3 + rl) instead of Eq.©. We try 
out this model in order to study the halo model de- 
pendence. For A1689 po changes by a factor ~ 5, 
and hence it depends on the halo profile. However 
for other clusters po changes slightly (less than a 
factor 2), because there are no data point near the 
core radius ro and only Mq is determined through 
the amplitude of the shear. 

In Fig|6]we show the density po vs. the temper- 
ature T to put a constraint on the neutrino mass. 
The dashed lines correspond to neutrino mass in 
Eq.©. From the figure, the minimum neutrino 
mass is 2 - 3 eV for CL0024 and CL1358. The 
above results are c onsistent with the previous X - 
ray measurements (panders 2003t ISanders 20071 ). 
Since the current limit is m„ < 2 eV from tritium 
/? decaj@, these values are comparable to or larger 
than this limit. 



"Sanders (2007) recently revised his calculation and gave ~ 
1/3 times smaller density than Sanders (2003). We confirm 
his calculation. Then the minimum neutrino mass is ~ 
3 1//4 = 1.3 times larger. Although Sanders (2007) is not 
yet published and there may be some ambiguities about 
the factor, we adopt his revised model. 

9 Particle Data Group Home Page : http://pdg.lbl.gov/ 
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The mean density inside the innermost, not sec- 
ond innermost, data point is much higher than p2- 
The neutrino masses are 5 — 6 eV for CL0024 and 
3—5 eV for CL1358 in this case. Hence there is 
an ambiguity about the definition of the central 
density for these clusters. 

The core density of A1689 is highest and the 
minimum neutrino mass reaches 4—6 eV. To check 
our result, we compare the core density with the 
previous studies in GR. Since ~ 10 x g for 
r < tq from Fig^b), GR is valid near the core. 
Halkola et al. (2006) gave central mass distribu- 
tion by analyzing 107 multiple images of 32 back- 
ground galaxies. Their mass distribution is con- 
sistent with Broadhurst et al. (2005b) (see Fig. 17 
of Halkola et al. (2006)). From a velocity disper- 
sion of 1450 km/s and the core radius 77 kpc, the 
core density is p c — 0.01M o /pc 3 . This is roughly 
consistent with our result. 

Allen (1998) suggests that the lensing core mass 
is generally a few times larger than X-ray core 
mass for non-cooling flow clusters (see also Clowe 
& Schneider 2001 in the case for A1689). Some 
ideas are proposed to explain the discrepancies 
: clusters are not in dynamical equilibrium, non- 
thermal pressure such as turbulent and magnetic 
pressure plays an important role, elongation of the 
cluster or substructures along a line-of-sight (e.g. 
Hattori, Kneib & Makino 1999). In fact, Lokas 
et al. (2006) show that A1689 is surrounded by 
a few substructures aligned along a line-of-sight, 
by studying the galaxy kinematics. CL0024 has 
a second mass clump which is separated at 3' 
from the center and has 30% of total cluster mass 



(jKneib et al. 20031 ). Jeeetal. (2007) recently sug- 
gest that CL0024 would be the merging cluster in 
line-of-sight direction. These systematics would 
affect our results and change the neutrino mass 
limit by factor (since m v is not very sensitive to 
Pc m„ oc p X J 4 from Eq.®). 

5. Results by Other Interpolation Func- 
tions 

So far, we used a standard interpolation 
function jl(x) — x/yl + x 2 alone. However 
the standard p is not consistent with TeVeS 
(jBekenstein 20041) . In this section, we also ex- 
amine other interpolation functions, Bekenstein's 
toy model in TeVeS p(x) = 4x/(l + y/1 + Ax) 2 
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Fig. 5. — The mean shear of 42 SDSS clusters. 
The solid and dotted lines are the beta models, 
and the dashed line is the dark halo model. 



10- 1 



o 

Q. 



10^ 



--"""~8eV 
_„---~"~6eV 

4eV 


A1689 
CL1358 


3eV , 






CL0024 


2eV_„ 





10 3 4 5 6789 10 

T(keV) 
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ter corresponding to the two density profile mod- 
els in Table. 1. The vertical dotted line connecting 
these pairs represents the systematic uncertainty 
in the density model. The dashed lines correspond 
to minimum neutrino mass in Eq.®. 
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Table 1 

Best fitting model for dark halo profile and central density 





M (M ) 


r (kpc) 


Po (M /pc 3 ) 




A1689 


(1.1 ±0.1) x 10 15 
(5.1 ±0.3) x 10 14 


174 ±11 

239 ±9 


(5.1 ±1.0) x 10' 
(9.0 ± 1.1) x 10" 


-2 
-3 


CL0024 


(2.0 ±0.4) x 10 14 
(1.7 ±0.2) x 10 14 


42 ± 15 
147 ±19 


(2.6 ±0.6) x 10- 
(2.8 ±0.3) x 10- 


-4 
-4 


CL1358 


(6.9 ±2.8) x 10 ia 
(4.0 ± 1.1) x 10 13 


69 ±31 
127 ±27 


(1.7 ± 1.0) x 10" 
(1.9 ±0.7) x 10" 


-3 
-3 



Note. — The best fitting parameters Mq and rn. The density 
po is the core density p c for A1689 and the mean density inside 
the second innermost data point p2 for the others. The upper 
row is the case of the halo profile M(< r) = M^r 3 /(r + rn) 3 
given in Eq.©, while the lower row is M(< r) = Mqt 3 /(r 3 + 

r 3 )- 



and simple model p(x) — x / (1 + x) in Famaey & 
Binncy (2005), in order to study the robustness of 
our results. Since the clusters have <?jv ~ 9o from 
Figs|2]|4] (b), our conclusions may depend on the 
choice of jl. 

In FiglJl we show the MOND predication for 
A1689 for the three types of p(x). The solid line is 
the standard model, the dashed line is the Beken- 
stein's toy model (B04), and the dotted line is the 
Famaey & Binney's model (FB). The standard 
model predicts the lowest value, because p(x) is 
highest at x ~ 1. B04 and FB show several times 
10% higher values than the standard model. 

The minimum neutrino mass only ~ 10% 
changes and the choice of p is not crucial. 

6. Conclusion 

We have studied the weak lensing of galaxy 
clusters in MOND. We calculate the shears and 
the convergences of the background galaxies for 
three clusters (A1689, CL0024, CL1358) and the 
mean profile of 42 SDSS clusters, and compare 
them with the observational data. It turns out 
that the MOND cannot explain the data irrespec- 
tive of go unless a dark matter halo is added. We 
also examine the three types of interpolation func- 
tion, but the conclusion docs not change. The 







— standard 






— B04 






FB 










. A1689 


4 





10° 10 1 

(arcmin) 



Fig. 7.— The MOND predications for A1689 for 
the three interpolation functions. The solid line 
is the standard model p{x) = x/yl + x 2 , the 
dashed line is the Bekenstein's toy model (B04) 
p{x) = 4x/(l + VI + 4x) 2 , and the dotted line 
is the Famaey & Binney's model (FB) p(x) = 
x/(l + x). 
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above results are consistent with those of previous 
studies (e.g. Aguirre, et al. 2001; Sanders 2003). 
If the dark halo is composed of massive neutrinos, 
its minimum mass is 4 — 6 eV for A1689 and 2 — 3 
cV for CL0024 and CL1358. However our results 
still depends on the dark halo model and inner 
data points. Even for A1689, the system with the 
most constraining data, the systematic is more im- 
portant than random errors as shown in Fig[6l In 
addition, there are some systematic uncertainties 
such as an elongation of cluster along line-of-sight. 
These effects would reduce the minimum mass by 
a factor of ~ 2. In the low acceleration region 
(gN < go), the external gravitational field becomes 
i mportant and would affect the outer shear p rofile 
(|Bekenstein fe Milgrom 1984 IWu et al. 20071) . 

In conclusion, the more careful study is neces- 
sary to put a stringent constraint, for example, a 
combination of the weak and strong lensing, X-ray 
gas and galaxy dynamics. However, it is beyond 
the scope of this paper and we will study as an fu- 
ture work. Even so, at present, we find that there 
is some tension between the lower bound of neu- 
trino mass in neutrino dark halo model in MOND 
and the upper bound by experiments. 

We would like to thank Takashi Hamana for 
useful comments and discussions. We also thank 
the anonymous referee for helpful comments to 
improve the manuscript. TC was supported in 
part by a Grant-in-Aid for Scientific Research 
(No.17204018) from the Japan Society for the Pro- 
motion of Science and in part by Nihon University. 
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